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Abstract We present the first simultaneous trapping 
of two different ultracold atomic species in a conserva- 
tive trap. Lithium and cesium atoms are stored in an 
optical dipole trap formed by the focus of a CO2 laser. 
Techniques for loading both species of atoms are dis- 
cussed and observations of elastic and inelastic collisions 
between the two species are presented. A model for sym- 
pathetic cooling of two species with strongly different 
mass in the presence of slow evaporation is developed. 
From the observed Cs-induced evaporation of Li atoms 
we estimate a cross section for cold elastic Li-Cs coUi- 
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Mixtures of ultracold atoms add many new possibilities 
to the study of interactions in ultracold gases since the 
different components of the mixture can have widely dif- 
ferent intrinsic properties, e.g., different masses, spins, 
scattering properties and quantum statistics. On one 
hand, many interesting phenomena are exclusive to gas 
mixtures, such as mixing and de-mixing phenomena ||l|, 
0], the subtle interplay between bosonic and fermionic 
atoms [y[ , or the prospect of forming ultracold heteronu- 
clear molecules through photoassociation Q . On the other 
hand, one may see an atomic gas mixture as a tool for 
transferring the favorable thermodynamic properties of 
one component (low temperature, high magnetization) 
to the other component. In this case one would like to 
use one gas that is easy to cool or polarize, and bring it in 
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contact with a second gas that is harder to handle. These 
processes are commonly referred to as sympathetic cool- 
ing and polarization exchange. 

Already 15 years ago, cryogenically cooled helium 
was used as a buffer gas to precool and load atomic 
hydrogen into a magnetic trap [Q. More recently, this 
method was extended to other atomic species and even 
molecules Q. Sympathetic cooling with ultracold, trapped 
gases has initially been studied for different hyperfine 
states of evaporatively cooled rubidium atoms in or- 
der to create two overlapping Bose-Einstein Condensates 
(BECs) m. Later, evaporative cooling of a mixture of 
magnetic spin states was employed as the key step to 
create a degenerate Fermi gas of fermionic '^'^K S . Most 
recently, sympathetic cooling in magnetic traps has been 
applied to different isotopes of the same atomic species: 
the bosonic species ^^Rb cooled by ^^Rb pi and the 
fermionic ^Li cooled by ^Li [llOllO]. So far, no experi- 
ments are reported which investigate elastic interactions 
in two distinct atomic species. 

Sympathetic cooling of a mixture of different atomic 
species in principle offers an enormous advantage if one 
species can be optically cooled to much lower tempera- 
tures than the other. We are investigating a mixture of 
the alkaline atoms lithium and cesium. This combination 
represents the prototype of a mixture in which one com- 
ponent (lithium) is hard to cool optically while the other 
(cesium) allows sub-Doppler temperatures ( < 3 /iK in 
free space) to be reached by simple polarization-gradient 
cooling. In contrast to evaporative cooling, the cooling 
process does not lead to a significant loss of either the 
cooled gas or the cooling agent. Sympathetic cooling 
with cesium (or other sub-Doppler cooled atoms) may 
be the method of choice for cooling many other atomic 
and molecular species, including LiCs molecules formed 
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by photoassociation from the same mixed-species sam- 
ple! 

We present the first observation of simultaneous trap- 
ping of two different atomic species in a conservative 
trap: lithium and cesium atoms stored in a very far- 
detuned optical dipole trap. Our optical trap is based 
on a CO2 laser emitting light at 10.6 /im [O]. The trap 
provides a conservative potential which is independent of 
the hyperfine state and spin projection of the atoms |18| . 
This allows us to study exoergic spin exchange collisions 
as well as elastic interactions. Exoergic collisions be- 
tween different species have been previously investigated 
in light-pressure traps |lJ,|l^;Q but the near-resonant 
light is known to influence the cross-sections even for 
collisions involving only the ground state |15[| . No previ- 
ous measurements of two-species collision cross-sections 
in conservative traps exist. If all atoms are in the lower 
hyperfine state, exoergic collisions are impossible. We 
nevertheless see slow evaporation of lithium, which is an 
indication that thermalization through elastic collisions 
is taking place. To interpret this evaporation and to es- 
timate a collison cross section from it, we adapt a model 
given in Ref. u& and we develop a rate equation model 
for sympathetic cooling with slow evaporation. 

This paper is organized as follows: An overview of 
the apparatus is given in section and the methods for 
loading and characterization of the trap are discussed 
in section pi The exoergic spin-changing collisions are 
studied in section m Finally, a rate equation model for 
sympathetic cooling is given in section ra, and compared 
to numerical results and experimental data. 



2 Quasi-electrostatic trap 

The focus of a CO2 laser beam constitutes an almost 
perfect realization of a conservative trapping potential 
for atoms. The atoms experience an optical dipole force 
which points towards the maximum of the intensity thus 
providing stable confinement in all spatial directions. 
Since the frequency of the CO2 laser is far below any 
dipole-allowed atomic transition frequency, the induced 
atomic dipole moment oscillates in phase with the os- 
cillating electric field of the laser as if the field was a 
slowly varying static field. Therefore, the name quasi- 
electrostatic trap (QUEST) has been coined for this par- 
ticular realization of an optical dipole trap. Due to the 
large detuning of the laser frequency from any atomic 
resonances, heating through photon scattering can be 
completely neglected. 

For an intensity distribution /(r) the potential C/(r) 
is given by C/(r) = -Q;stat^(r)/2£oc where ofstat denotes 
the static polarizability of the atoms |l^,ttq]. For a fo- 
cused beam of power P and waist w one gets a trap depth 



et = 2asta,tP/ {tteocw'^) . We employ an industrial, sealed- 
tube CO2 laser (Synrad Evolution) providing 140 W of 
power in a nearly TEMqo transversal mode character- 
ized by Af^ = 1.2. The laser beam is first expanded by 
a telescope, and then focused into the vacuum chamber 
by a lens of 254 mm focal length. The focus has a waist 
of 85 /im, yielding a trap depth of 390 /iK for Li and 
1000 /iK for Cs (in units of the Boltzmann constant fee). 
Assuming a pure Gaussian beam, one finds a Rayleigh 
range zr of 1.5 mm. The axial and radial oscillation fre- 
quencies in the harmonic approximation are given by 
Wz = (2C/o/™-Zr)^^^ and lo,- = [AUo/mw'^Y/'^ with m 
denoting the mass of the atom. 

Fig. P shows a schematic view of the apparatus. Atoms 
are transferred into the dipole trap from magneto-optical 
traps (MOT) for lithium and cesium which are superim- 
posed on the focus of the CO2 laser beam. Both MOTs 
are loaded from atomic beams that are cooled and decel- 
erated by decreasing-field Zeeman slowers which use the 
fringe fields of the MOT magnetic quadrupole field for 
the final deceleration stage. Fluorescence of the atoms 
in the MOT is detected by photodiodes with appropri- 
ate filters to discriminate the resonance fluorescence of 
Li at 671 nm from the fluorescence of cesium at 852 nm. 
The main vacuum chamber at a background pressure of 
about 6 X 10~^^ mbar is connected to the oven chambers 
by tubes which are divided into differentially-pumped 
sections. Both atomic beams can be interrupted by me- 
chanical shutters. 



3 Storage in the quasi-electrostatic trap 

3.1 Cesium 

To load cesium into the QUEST, a cloud of typically 10^ 
atoms is flrst accumulated in the MOT with a quadrupole 
fleld gradient of 26 G/cm along the symmetry axis. The 
atoms are further compressed and cooled by decreasing 
the detuning of the MOT trapping laser from initially 
-2rcs to -25rcs with respect to the Q'^Si/2{F = 4) - 
6^P3/2(-P' — 5) transition (natural linewidth /cs/27r — 
5.3 MHz). After 5 ms of compression, the distribution of 
atoms in the MOT has a rms radius of 400 /im corre- 
sponding to an mean density of about lO^*' atoms/cm^. 
After the compression, the magnetic fleld of the MOT 
is turned off, and the atoms are further cooled in a blue 
optical molasses ]1^ by detuning the laser field to -l-2/cs 
with respect to the 6^81/2 (F = 4) - 6^P3/2(-F = 4) tran- 
sition for 20 ms. 

The CO2 laser beam is present during the whole load- 
ing phase. Since the static polarizability of the 6^P3/2 
state in cesium is larger than the polarizability of the 
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Fig. 1 Schematic setup of the experiment. The insets show the relevant energy levels of lithium and cesium. 



ground state by around a factor of 4, the atomic res- 
onance frequency of the cesium D2 line is inhomoge- 
neously shifted with a maximum value of roughly — 10/cs- 
However, the number of atoms in the MOT is found to 
be independent of whether the CO2 is turned on or off. 

After the cooling laser beams have been turned off, 
up to 10^ atoms remain trapped in the focus of the CO2 
laser beam at a temperature of 30 /iK and a peak density 
close to 10"'^^ atoms/cm'^. Atoms are prepared in either 
the F = 3 or the F = A hyperfine ground state by block- 
ing the cooling laser 1 ms after or before the repumping 
laser has been blocked, respectively. To detect the op- 
tically trapped atoms, we either recapture the atoms in 
the MOT and measure the intensity of the fluorescence 
light, or we analyze the spatially resolved absorption of a 
weak resonant probe beam passing through the trapped 
atoms (absorption imaging |£^). The latter method al- 
lows one to absolutely calibrate the atom number. In 
Fig. pi the absorption image of the trapped Cs atoms 
shows the expected elongated shape with an aspect ra- 
tio of tJax/wrad — 44. Notc, that the radial extension of 
the atomic cloud is not resolved by our camera system. 

The evolution of the number of atoms in the dipole 
trap is measured by recapturing the atoms back into the 
MOT after a variable storage time. The number of atoms 
is determined from a charge coupled device camera pic- 
ture by integrating over the fluorescence image, while 




Fig. 2 Absorption image of cesium atoms in the QUEST. 
The dashed lines indicate the profile of the CO2 laser beam 



subtracting a background image. We expect the parti- 
cle number determined by this method to be accurate 
within a factor of two. Fluctuations in the loading flux 
are cancelled by loading the MOT up to a certain density 
(determined by the fluorescence signal) rather than for 
a fixed period of time: A photodiode records the fluores- 
cence of the atoms in the MOT and triggers the transfer 
cycle as soon as it exceeds a predefined discriminator 
level. The decay of the number of stored cesium atoms 
is shown in Fig. H(a) from which one infers a 1/e storage 
time of 150 s. The very long storage time is purely lim- 
ited by collisions with hot atoms out of the background 
gas at the base pressure of6x 10~^^ mbar in our vacuum 
chamber p3|. Laser-noise induced loss mechanisms [El] 
can apparently be completely neglected. Since Cs has a 
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Fig. 3 Very long storage times for (a) cesium and (b) 
lithium in the quasi-electrostatic trap. The atoms are opti- 
cally pumped into the lowest ground state. The correspond- 
ing 1/e lifetimes are 150 ±15 s for cesium and 125 ± 15 s for 
lithium. 



large scattering length the Cs gas thermalizes on a time 
scale of < 0.1 s. In contrast to previous experiments in 
a shallower trap [ p2| , spontaneous evaporation of atoms 
is exponentially suppressed by a factor exp(— 77), where 
ry = et/ksT is more than 20. Experimentally, we indeed 
observe neither loss of Cs atoms nor reduction of the 
temperature due to evaporation. 

To determine the characteristic frequencies of the 
trap, we induce axial and radial oscillations of the cloud 
of trapped atoms. Radial oscillations are excited by a 
release-recapture cycle in the optical trap. First, the 
laser is turned off (switch-off time < 100 /is) which leads 
to a ballistic expansion of the atomic cloud. After ^^ 
0.5 ms, the cloud has expanded to about twice its initial 
radial width. Then, the CO2 laser is turned on again 
and almost all atoms are recaptured by the dipole trap. 
Through this process, the atoms have gained potential 
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Fig. 4 Measurement of (a) radial, and (b) axial oscillations 
of Cs atoms in the quasi-electrostatic trap. The solid lines 
show the fits to a damped sinusoidal oscillation. For the axial 
oscillation, a constant drift of the focus position was included 
into the fit. 



energy, which is subsequently converted into kinetic en- 
ergy. Thus, the radial size of the atom cloud starts to 
oscillate at twice the radial oscillation frequency. Since 
we can not fully resolve the radial extension of the atom 
cloud with our absorption imaging system, we measure 
the corresponding oscillation in the kinetic energy. For 
this purpose, the CO2 laser is turned off again after 
a variable delay relative to the release-recapture cycle, 
and the radial extension is determined through absorp- 
tion imaging after 1.2 ms of ballistic expansion. Fig. |4|(a) 
shows the oscillation of the root mean square velocity at 
a frequency of 1.6 kHz. The corresponding oscillation fre- 
quency of the trapped atoms of 0.8 kHz is in reasonable 
agreement with the expected radial oscillation frequency 
for a focused beam with a waist of 85 /im. 

Axial oscillations can easily be excited by loading at 
an axial position off the focus. For this purpose, the rela- 
tive position of the MOT is shifted along the propagation 
direction of the C02-laser beam. After transfer into the 
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trap depth radial frequency axial frequency 



Lithium 390 fiK 50 Hz 

Cesium 1000 ^lK 18 Hz 



2.2 kHz 
0.8 kHz 



Table 1 Parameters of the quasi-electrostatic trap for 
lithium and cesium. The Li frequencies are calculated from 
the measured Cs frequencies. 



QUEST, the cesium cloud starts to perform a damped 
oscillation at a frequency of 18 Hz as depicted in the up- 
per graph of Fig. ^ The damping of the oscillation is 
caused by the anharmonicity of the confining potential. 
Assuming a pure Gaussian beam, one would expect an 
oscillation frequency about twice as large as the actual 
measured value. We explain the discrepancy by aberra- 
tions in the optical system. Table |3.l| summarizes the 
relevant trap parameters for lithium and cesium. 

In the first few seconds after the CO2 laser has been 
turned on, we observe a small drift of the axial position 
of the laser focus, which can be attributed to thermal 
effects in the laser cavity and optics. To eliminate this 
drift of the focal position, we leave the CO2 laser perma- 
nently on while loading the MOT and transferring the 
atoms. 



3.2 Lithium 

For cesium the thermal energy of the atoms in the MOT 
is much smaller than the depth of the QUEST which re- 
sults in a rather large transfer efficiency. The situation 
is different for lithium. Due to the large recoil of pho- 
tons on lithium atoms, the thermal energies in a lithium 
MOT are comparable or even larger than the depth of 
the dipole trap resulting in a poor phase-space overlap 
and thus a small transfer efficiency. 

The lithium MOT is also loaded from a Zeeman- 
slowed atomic beam. Details of the MOT system are de- 
scribed in Ref. ||2^ . The largest numbers of Li atoms are 
accumulated at relatively large detunings of the MOT 
lasers acting on the two hyperfine components of the Li 
D2 line. After loading typically lO'' atoms into the MOT, 
the atomic cloud is compressed for 4 ms by switching 
the detunings closer to resonance (—1 Fu with Fu/^tt — 
5.9 MHz). Due to the unresolved hyperfine structure of 
the excited state, sub-Doppler cooling mechanisms are 
inefficient for lithium and temperature is essentially de- 
termined by the limit for Doppler cooling (140 ^K at low 
saturation). To attain a temperature close to that limit, 
the MOT magnetic field is turned off and the MOT lasers 
are attenuated to 20% of their full power which leads to 
a final temperature of 350 fiK and a cloud size of 800 /xm 
after typically 1 ms. In contrast to cesium, the light shift 



of the CO2 laser beam on the D2 resonance line is neg- 
ligible since the static polarizability of the ground and 
excited state differ by only 25%. 

After all MOT laser beams have been shut off, about 
0.5% of the Li atoms are captured by the QUEST with 
a peak density of ~ 10^° atoms/cm"^. Thermalization of 
these atoms is extremely slow: estimated collision times 
based on the s-wave scattering lengths for Li are exceed- 
ing 100 seconds. On the time scales of the experiments 
discussed in the following, lithium can therefore not be 
assumed to be in thermal equilibrium. 

Due to the small optical density of the sample, we 
were so far not able to directly measure the distribution 
of lithium atoms in the trap by absorption imaging tech- 
niques. Instead, we have performed a Monte-Carlo simu- 
lation of the transfer process which provides us with the 
spatial and energy distributions. The energy distribu- 
tion is highly non-thermal with most atoms populating 
high lying energy states in the potential because of the 
large initial thermal energy relative to the trap depth. 
The small rate for elastic collisions inhibits evaporation 
of these loosely bound atoms. The measurement of the 
number of trapped lithium atoms as a function of stor- 
age time depicted in Fig. ^ thus shows a pure exponential 
decay with a time constant of 125 s with no indication 
for evaporation of particles from the trap. On the ba- 
sis of the above estimate for the thermalization times, 
we assume the energy distribution hardly changes with 
time. 



3.3 Combined transfer 

To load a mixture of lithium and cesium, one could 
start from a two-species MOT jlj], but in this case the 
numbers of atoms and the densities are severely lim- 
ited by light-assisted inelastic collisions between the two 
species. However, our measurements with a 2-species 
MOT showed that collisions involving optically excited 
Li (Li*) and ground-state Cs are strongly suppressed by 
optical shielding through the repulsive van-der-Waals in- 
teraction in the Li*-Cs manifold |jlj|. It is therefore fa- 
vorable to first optically cool and transfer Cs into the 
QUEST, and then load the Li MOT. A typical loading 
time of the Li MOT is 4 s. In order to exclude any loss 
processes by inelastic collisions, we spatially separate the 
Li MOT from the optically trapped Cs sample with a 
small bias magnetic field (3 G) shifting the centre of the 
Li MOT by about 2 mm. In this way, we observe no loss 
of Li and only a 50 % loss of Cs in the combined transfer, 
as compared to the separate trapping of only one species. 
A detailed investigation of the combined transfer will be 
published elsewhere. 
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4 Spin-changing lithium-cesium colHsions 

Two atoms approaching each other in their electronic 
ground state may either colhde elasticaUy with the to- 
tal energy being conserved, as will be discussed in the 
next Section, or they may undergo an exoergic collision 
transforming internal energy into kinetic energy of the 
collision partners. The latter process leads to trap loss 
if the kinetic energy release is larger than the escape en- 
ergy of the trap. Optical dipole traps are perfectly suited 
to study inelastic processes among trapped atoms since 
the optical trapping force is conservative with a well- 
defined potential, the investigation may be performed in 
zero magnetic field or in a controlled homogeneous field, 
and the atoms can be prepared in any internal state by 
appropriate optical pumping methods. Up to now, no 
experiment on ultracold inelastic collisions between two 
different atomic species in the electronic ground state 
has reported. 

For ground-state Li and Cs atoms in an optical dipole 
trap, the essential inelastic processes are collisions that 
change the hyperfine quantum numbers F and mp . The 
spin-changing collision process between lithium and ce- 
sium can be written symbolically as 

Li(F=l) + Cs(F = 4) -> UCsCT^g or ^S^) 

^ Li(F=l) + Cs(F = 3) + ziekin • 

The potential curves for the singlet and triplet ground 
states of the Li-Cs pair are schematically depicted in 
Fig. y. Spin-changing collisions exiting at a lower en- 
ergetic asymptote than the entrance channel lead to a 
transformation of the internal energy into kinetic energy 
Z\ekin of the collision partners, and consequently to loss 
from the trap. Although lithium takes the major share 
of the kinetic energy gain, both atoms escape from the 
QUEST due to the small trap depth as compared to the 
hyperfine energy. 

Inelastic collisions of trapped ground-state atoms have 
been extensively studied theoretically and experimen- 
tally since the early days of magnetic traps [|4| . Hyperfine- 
changing collisions represent the major limitation on achiev- 
able densities and/or storage times in magnetic traps. 
These collisions in magnetic traps are almost unavoid- 
able, since atoms can not be trapped in the lowest ener- 
getic state. In particular, spin-changing collisions impede 
sympathetic cooling between different atomic species in 
magnetic traps. 

Exoergic ground-state collisions also turned out to be 
an important loss mechanism in shallow light-pressure 
traps. They were observed and analyzed for homonu- 
clear collisions p^ , as well as for mixtures of different 
elements Jl3,0,ll3. However, as has been pointed out 
recently [|15[, the presence of near-resonant light strongly 




Fig. 5 Ground-state potential energy curves for the het- 
eronuclear dimer LiCs. The inset shows an enlargement at 
long range. The separated atom hyperfine levels are char- 
acterized by the total angular momentum F {Fcsi= 3 and 
Fcst= 4 for Cs, -FLi|= 1 and Fuf= 2 for Li). The arrow 
indicates the incoming channel. 



modifies the dynamics of the collision process which can 
lead to a dramatic change in the measured rate coeffi- 
cients. Therefore, care should be taken to interpret trap 
loss and the corresponding rate coefficients in shallow 
light-pressure traps as an indication of pure ground- 
state collisions. In a purely conservative trap such as the 
QUEST the situation is much better defined since the 
trapping light has negligible influence on the collision 
process. 

To qualitatively investigate whether hyperfine-changing 
collisions between Li and Cs can be observed with the 
particle numbers and densities available in our QUEST, 
we have pumped one of the species into the upper hyper- 
fine ground state and recorded the effect on the storage 
time of the other species, which was prepared in the 
lower hyperfine state. The result of these experiments is 
depicted in Fig. |6[ From the timescale of the graph it 
immediately becomes apparent, that the presence of a 
hyperfine-excited atomic gas reduces the storage time of 
the species in the lowest energetic state by more than an 
order of magnitude as compared to the rest-gas limited 
lifetimes (see Fig. ||). 

In Fig. 0(a), we show data for a combination of Li in 
the upper hyperfine state [Li 2^5i/2(-FLij,=2)], with Cs 
in the lower hyperfine ground state [Cs 6^51/2 (-FcsT^ 
3)]. Without the presence of Cs, the Li slowly decays 
through exoergic Li-Li collisions. The effect of collision 
with restgas atoms can be neglected on this timescale. 
In presence of Cs, the Li decay is much faster. Spin- 
exchange between Li and Cs atoms is the dominant pro- 
cess for these interspecies collisions. An enhanced Li-Li 
inelastic collision rate due to an increased Li density af- 
ter thermalization with the Cs may also play a minor 
role, but this would lead to a non-exponential decay of 
the Li, while we observe an exponential decay with a 
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Fig. 6 Exoergic spin-changing collisions between lithium 
and cesium resulting in trap loss, (a): Li 2^51/2 (-Flit =2) 
and Cs 62S'i/2(Fcsi=3). (b): Li 2'^S',j2{Fui^l) and Cs 
6^>S'i/2(fcsT=4). Here, | stand for the upper hyperfine 
ground state and [ stands for the lowest-energy hyperfine 
state. The dashed lines indicate the fit to an exponential de- 
cay yielding a time constant of 1.7 s and 1.9 s for the upper 
and lower graph, respectively. 



time constant of 1.7 s. The effect of spin-changing Li-Cs 
collisions on the Cs is not observable as the Cs atoms 
outnumber the Li atoms. As the lithium distribution is 
not known, and may even be nonconstant due to elas- 
tic collisions happening simultaneously, we can currently 
only estimate the order of magnitude of the inelastic rate 
constant /3csiLiT '^ 10~^° cm^/s. We are currently im- 
plementing an improved absorption imaging system for 
Li in our apparatus, which will allow us to measure the 
density distribution of the Li atoms und thus allow for 
quantitative measurements of the rate coefficients. 

In Fig. P(b), we investigate the influence of Cs in 
the upper hyperfine ground state [Cs 6^S'i/2(-F'csT=4)] 
on Li in the lower hyperfine state [Li 2^S'i/2(-FLix=l)]. 
Cs itself decays due to exoergic Cs-Cs collisions. The Li, 



however, is lost on a much shorter timescale of 2 s. Due 
to the unknown Li distribution, we can only offer an 
estimate of the rate constant of the exoergic collisions 
of /JcsTLii =^ 10^^" cw? /s. As one might expect, both 
rate constants are of the same order of magnitude. 



5 Sympathetic cooling with slow sympathetic 
evaporation 

In case both atomic gases are in the lower hyperfine 
states [Li 2'^Si,2{Fui = l) and Cs 62S'i/2(Fcsi=3)], only 
elastic collisions are possible. These collisions lead to en- 
ergy transfer between the distributions, which in turn 
leads to two different phenomena: The most obvious 
one is the equalization of the temperatures of the two 
gases, which is referred to as sympathetic cooling. A sec- 
ond phenomenon is the evaporative loss of one or both 
species through collisions with the other species, which 
we call sympathetic evaporation. 



5.1 Sympathetic cooling 

For the case where all relevant temperatures are well 
below the trap depth, the dynamics are relatively sim- 
ple: No atoms evaporate and the temperatures of the two 
components approach each other. A quantum-mechanical 
treatment of sympathetic cooling has been given in Refs. [|26 
p7| . For the Boltzmann regime, an analytical formula for 
the evolution of the temperatures in a mixture of equal 
mass atoms was derived in Ref. |lq]. In the same paper, 
the approach is extended towards atoms with different 
mass, however, the dependence of the cooling efficiency 
on the mass ratio is derived under the assumption of 
equal trap frequencies. This assumption corresponds to 
an implicit mass dependence of the trap potential, which 
is not present for realistic traps, and which leads to a 
strongly overestimated reduction of the cooling efficiency 
at large mass ratio. Below, wc adapt the approach of 
Ref. |lq] by assuming an explicit trap potential, so that 
our results remain valid for large mass ratios. The aver- 
age energy transfer per collision can be calculated from 
simple kinematics. 



AE^^2^ikB{Ti-T2); ^ 



Amim2 
(mi + 777,2)^ 



(1) 



The energy transfer is reduced by a factor ^ with re- 
spect to collisions between atoms with equal mass. For 
equal mass atoms in a harmonic trap, approximately 3 
collisions are needed for thermalization (see [M and ref- 
erences therein), as the heat capacity per atom is Sfc^. 
For ^■^^Cs and ^Li, we need 3/^ « 15 collisions to ther- 
malize. 
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The Cs-Li collision rate depends on the spatial over- 
lap between the thermal distributions, 



rcoii = 0-12V / (i^xni(x)n2(x) 



(2) 



Here v is the mean thermal relative velocity, 



-,1/2 



TT mi 7712 



(3) 



We assume a harmonic trapping potential, 14 (x) = (a,; • 
x)-^, where i = 1, 2, and cti is a species-dependent vector 
of force constants. We will assume the trap constants are 
linearly dependent, i.e. a2^Pai=f3a. This allows us to 
write the collision rate as 



rcoii — 



(T12VN1N2 



nUyaz (Ti + T'2/3 



-2\-3/2 



(4) 



The timescale of thermalization can be obtained from 
the differential equation for AT = T2 — Ti : 



_ _i _ djAT) _ 4(JVi + N2)cTi2S.a^ayaz ^, ,^. 



ATdt 



SiT^k, 



(Ti/2mi+r2/2m2)i/' 

(ri + r2/3-2)3/2 

For equal masses and trap frequencies we reproduce the 
resuhof Ref. |§. 

For the initial conditions of our experiment we obtain 
an approximate initial cooling time constant, where we 
neglect the influence of the Cs temperature since Cs is 
much colder, heavier and more tightly confined than the 
Li. 



"' " 37r2fcBrLi(2mLi)i/2 



(6) 



In the final stage of sympathetic cooling the final tem- 
perature T = (iViTi + N2T2)/{Ni + N2) is approached 
exponentially, with a time constant 



4V2g(A^i + N2)cji2axayaz 
2,'n'^kgf{mi+m2Y''^ 



(l+^-2)-3/2_(7) 



The effect of the mass difference on the thermalization 
time is a factor ^"^/^ ~ 2.3 in the time constant. At 
large mass difference, the thermalization efficiency per 
collision is reduced by a factor ^, but the collision fre- 
quency depends on u, which is proportional to ^^^'"^ . 



5.2 Sympathetic evaporation 

Under the circumstances in our experiment, we expect 
evaporation of the Li atoms especially during the first 
few collision times, as the initial Li energy distribution 
has a width comparable to the trap depth. Once the 
gases are thermalized, Li evaporation is exponentially 
suppressed, but still may occur due to the high Li-Cs 
collision rate. Cs evaporation is always strongly sup- 
pressed as the Cs trap depth is much larger than the 
relevant temperatures. We can approximately describe 
the regime of slow evaporation by assuming a truncated 
Boltzmann distribution pq ] for the Li atoms. The num- 
ber of evaporation events per second is the number of col- 
lisions multiplied by the probability Povap that a collision 
leads to evaporation. This probability can be found from 
a detailed-balance approach in the following way: If the 
gas is in full thermal equilibrium (i.e. at the same tem- 
perature and chemical potential) with a gas of identical 
untrapped atoms, evaporation will be exactly compen- 
sated by "atom capture" events. The untrapped atoms 
will have a density n„(r) which is to good approxima- 
tion n(r)77^/2 exp(— 77). In case rjS^ '3> 1, practically every 
collision of a trapped with an untrapped atom will lead 
to capture of the untrapped atom, i.e., the number of 
capture events is well approximated as the number of 
collisions between untrapped atoms and trapped atoms, 
which is a fraction rj exp(— 77) of the total number of colli- 
sions. It follows from the detailed balance approach that 
this equals the number of evaporation events, so that we 
find the following approximation: 



-Povap ~ 77exp(-77). 



(8) 



This evaporation probability leads to an evaporation 
rate of r'ovap = Aoi-Pcvap- The amount of energy lost 
in an evaporation event is of order {r]+ VjksT, while the 
Cs atom that is involved in the last collision only loses a 
small amount of energy, of order fc^T. The Li distribu- 
tion on the other hand loses a particle and an energy of 
r^ rjksT, leading to a decrease of the mean energy per 
Li atom of 

^(eLi) ^ iv - 3/2 - S) 

keT A^Li ■ ^ ' 

Here 5 represents the average potential energy, 5 — Q 
for a box potential and 5 = 3/2 for a harmonic trap. 
Since the evaporation cools the Li more than the Cs, 
we expect that evaporation leads to a small difference in 
mean energy between the two components, of order 



ecs 



eLi 



2-e 



'"(^-3/2-5-770. 



(10) 



Heree = (Afcsecs+^^LieLi)/(A'^Cs+^Li)- The evaporation- 
induced energy difference for a mixture of equal mass 
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Fig. 7 Evolution of the particle number and internal energy 
in a simulated thermalization experiment. The number of Cs 
atoms is constant at 5 times the initial number of Li atoms. 
The mean energies cross over after 45 collisions, see text. 
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Fig. 8 Simulated evolution of the Li and Cs distributions. 
The left pane shows the initial Li distribution, which is a 
strongly truncated Maxwell-Boltzmann distribution. The Cs 
distribution is shown reduced in height by a factor 5. 



atoms is of the same order of magnitude as the differ- 
ence in mean energy between a truncated Boltzmann 
distribution and the full thermal distribution described 
by the same temperature [ESl. For atoms with differ- 
ent masses, the thermalization is slower by a factor ^, 
whereas the evaporation rate is practically independent 
of the mass difference. This leads to an enhancement of 
the temperature difference in evaporating mixtures, as 
we have observed in simulated evaporation experiments. 
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Fig. 9 Evolution of the number of trapped Li (F=l) and Cs 
(F=3) atoms in simultaneous trapping experiments (closed 
symbols). For comparison, the open symbols show the evolu- 
tion of separately trapped gases. Thin lines: Exponential fits. 
Time constants are 150 s for Cs with or without Li, 125 s for 
Li without Cs, and 35 s for Li trapped simultaneously with 
Cs. 



in an analytical model. To obtain understanding of all 
phases of the thermalization process we made a numer- 
ical simulation of Li-Cs thermalization in a finite depth 
"box" potential well. We assume that due to the large 
scattering length the Cs distribution remains thermal at 
all times, while we numerically integrate the Boltzmann 
equation for the Li distribution. In Fig. |7| we show the 
evolution of the mean energies and particle numbers, 
where a broad initial distribution for the Li atoms was 
chosen. The evolution of the distributions is shown in 
Fig. |. 

In the initial phase of the thermalization process, the 
first ~ 5 collisions in Fig. |7|, the Li has an internal en- 
ergy per particle comparable to the trap depth, so that 
evaporation is very fast. After this rapid evaporation, a 
thermalization phase sets in (5 to 50 collisions) in which 
the mean energies rapidly converge. Finally, after 50 col- 
lisions a quasi-stationary phase with slow sympathetic 
evaporation and only slowly varying, almost equal tem- 
peratures sets in. We observe that although we start with 
Tu > Tcs, after ^ 50 collisions the evaporating Li gas 
has a lower average energy per particle than the Cs, rea- 
sonably well described by Eq.(nG). 



5.3 Numerical simulation 



5.4 Comparison with measurements 



The initial phase of thermalization, in which the distri- 
butions arc strongly non-thermal, is difficult to describe 



A complete characterization of the thermalization pro- 
cess requires measurement of all relevant thermodynamic 
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variables: the Li and Cs temperatures and particle num- 
bers. Unfortunately, presently the Li temperature cannot 
be measured, and as A^cs ^ ^Li, the effect of sympa- 
thetic cooling on the Cs temperature is too small to be 
measured accurately. 

The evolution of the number of Li atoms in the trap 
is shown in Fig. H. When stored separately, both Li and 
Cs gases decay due to rest gas collisions in 125 and 150 
s respectively. There is no indication of evaporation in 
either case: Cs does not evaporate as its temperature is 
far below the trap depth. Li evaporation is energetically 
possible but the cross section for Li-Li collisions is so 
small that the characteristic timescale is > 1000 s. 

When the gases are stored simultaneously, the Li 
atom number shows a non-exponential initial behavior 
followed by an approximately exponential decay with a 
characteristic time of 35 s, which is much faster than the 
125 s rest gas induced decay. We interpret the initial loss 
as due to evaporation in the first few Li-Cs collisions (cf. 
Fig. 1^), and the slow decay as due to sympathetic evap- 
oration after thermalization. Equation (0) now allows us 
to estimate the collision frequency. From the calculated 
Li trap depth and the measured Cs temperature after 20 
s, Tcs = 44(5) fjK we find Povap - 1-5 x IQ-^. This leads 
to an estimated collision rate per Li atom of ^ 20 s^^, 
and finally through Eq.M) we find the estimated cross 
section 

acsLi - 5 X lO-i^cm^. (11) 

From the collision rate we infer a cross-species ther- 
malization time of approximately 1 s, which is much 
shorter than the rest gas limited storage times in our 
trap. Therefore, we conclude that sympathetic cooling 
of Li by Cs is a viable method to produce ultracold Li 
gas. 



6 Conclusion 

In conclusion, we have simultaneously stored ultracold 
lithium and cesium atoms in a conservative optical dipole 
trap. Exoergic spin-changing collisions lead to trap loss 
on a time scale of a few seconds at densities of 10^" 
atoms/cm"^ and 10^^ atoms/cm"^ for lithium and cesium, 
respectively. Under conditions where exoergic collisions 
are impossible, a slow loss of Li atoms is observed, which 
we attribute to sympathetic evaporation through elastic 
Li-Cs collisions. The interpretation in terms of a model 
for sympathetic cooling and evaporation provides an es- 
timated cross-thermalization time of one second. 

The interpretation of the observations in terms of 
sympathetic cooling should be confirmed by measure- 
ments of the temperature evolution and phase space 
density. Such measurements are currently in prepara- 
tion. Temperature data will also enable us to measure 



the elastic cross section acsU for cold, elastic Li-Cs col- 
lisions more accurately. Optimization of the Cs cooling 
process is certainly possible and will lead to lower final 
temperatures of the Li gas. 

Repeated pulsed optical cooling of the Cs may be 
useful to remove thermal energy from the system with- 
out loss of atoms. This method offers a possible route 
towards quantum degeneracy of Li without evaporation. 
The de Broglie wavelength of Li is larger than that of 
Cs at the same temperature by a factor of 4.4, which 
makes it possible to achieve quantum degeneracy with 
either the fermionic or the bosonic Li isotopes in ther- 
mal contact with optically cooled Cs. For lO*' Li atoms 
in our trap we find a degeneracy temperature Tdog — 
h,{ujx(^y(^z ^LiY /ks = 3/zK, well within the range of 
standard optical cooling techniques for Cs. But even at 
the temperatures we reach presently, sympathetic cool- 
ing significantly increases the Li density in the centre 
of the trap. This is an essential step towards creation 
of cold hcteronuclear molecules, as all proposed forma- 
tion mechanisms, in particular photoassociation, rely on 
a high pair densities of the two species. 
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